Quantum limit of the triplet proximity effect in half-metal - superconductor junctions 
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We apply the scattering matrix approach to the triplet proximity effect in superconductor-half 
metal structures. We find that for junctions that do not mix different orbital modes, the zero 
bias Andreev conductance vanishes, while the zero bias Josephson current is nonzero. We illus- 
trate this finding on a ballistic half-metal-superconductor (HS) and superconductor-half-metal- 
superconductor (SHS) junction with translation invariance along the interfaces, and on HS and SIfS 
systems where transport through the half-metallic region takes place through a single conducting 
channel. Our calculations for these physically single mode setups - single mode point contacts and 
chaotic quantum dots with single mode contacts - illustrate the main strength of the scattering 
matrix approach: it allows for studying systems in the quantum mechanical limit, which is inacces- 
sible for quasiclassical Green's function methods, the main theoretical tool in previous works on the 
triplet proximity effect. 

PACS numbers: 74.50.-|-r, 74.45. +c,74.78.Na 



I. INTRODUCTION 

The recent experimental observation of the Joseph- 
son effect in a half-metallic junction between two 
superconducting reservoirsi has renewed interest in 
superconductor-ferromagnet hybrid devices. The obser- 
vation of a supercurrent in a half metal is remarkable, be- 
cause Cooper pairs in spin-singlet superconductors con- 
sist of a pair of electrons with opposite spin, whereas 
a half metal conducts electrons of one spin direction 
only^i^"^ The resolution of this apparent paradox is the 
so-called "triplet proximity effect", first predicted theo- 
retically by Bergeret, Volkov, and Efetov.^ (See also Ref. 

as well as Ref. for a review.) The triplet proximity 
effect relies on the conversion of spin-singlet Cooper pairs 
of electrons with opposite spin into pairs of electrons of 
equal spin at a spin-active interface between the super- 
conductor and the half metal.-'-!§ Since pairs of equal- 
spin electrons can be transmitted coherently through a 
half metal, the triplet proximity effect can indeed explain 
the observation of a Josephson current in the experiment. 

Most theoretical studies of the triplet proximity ef- 
fect were done using the quasiclassical Green's function 
methodi^i^i^iiSiiiii^iHaiiis rpj^jg method is appropriate 

for systems in which transport takes place through many 
conducting channels J^iii For systems with few channels 
only, the Green's function technique should be applied 
without the quasiclassical approximation. This, albeit 
doablcfi^ii^iii^ can lead to calculations of significant com- 
plexity. Another method that is particularly well suited 
for few channel structures is the scattering matrix ap- 
proach. This method has been frequently used in the 
context of transport problems involving superconductors 
(for a review, see Ref. [l^. However, it has not yet been 
applied to the triplet proximity effect. It is the goal of 
the present article to fill this gap. 

In the language of the scattering approach, the triplet 
proximity effect relies on the coherent Andreev reflec- 



tion of electron-like excitations into hole-like excitations 
with the same spin.^^ Conventional Andreev reflection, 
as it takes place at the interface between a normal metal 
and a superconductor, consists of the reflection of an 
electron into a hole with opposite spin. "Same spin" 
and "opposite spin" here refers to the spin band from 
which the electron and hole are taken. Since electron 
and hole from the same spin band have opposite angu- 
lar momentum, conservation of angular momentum im- 
plies that electron and hole are from opposite spin bands. 
Hence, Andreev reflection of electrons into holes from the 
same spin band requires that the interface between the 
half metal and the superconductor is spin active. Ex- 
amples of appropriate spin active interfaces are a thin 
ferromagnetic or half-metallic layer with a polarization 
that is non-coUinear with the half metal's polarization or 
a normal-metal spacer layer with strong spin-orbit scat- 
tering. 

Our focus is on systems with the fewest number of 
channels possible, a single conducting channel at the 
Fermi level. This limit can be achieved by having sin- 
gle channel contacts between the superconductor(s) and 
the half metal. As an example of this limit, we use the 
scattering theory to address the simplest single chan- 
nel half-metal-superconductor (HS) junction that can 
display triplet proximity effect: a single channel ferro- 
magnetic or half-metallic ballistic point contact between 
H and S electrodes. To study a more complex situa- 
tion, we investigate HS and superconductor-half-metal- 
superconductor (SHS) junctions where the half metal is a 
chaotic quantum dot with single-channel point contacts. 
We also study the case of ballistic devices which have 
translation invariance along the interfaces. This situa- 
tion allows for a single channel description as well, since 
the translation symmetry ensures that different trans- 
verse modes do not mix. While the latter system can in 
principle be addressed by the quasiclassical Green's func- 
tion method, the former, physically single channel setups 





t 






H 






s 




^—f' 







J 



D R 




R' C 



D R 



Figure 1: (Color online) Composite HS junction consisting 
of a half- metallic contact (left), a superconducting contact 
(right), and a spin-active intermediate layer (center). In most 
of our considerations, the intermediate layer is taken to be fer- 
romagnetic with a magnetization direction not coUinear with 
the polarization of the half metal. Transport through the 
HS junction is described by the scattering matrix TZ, which 
is calculated in terms of the Andreev reflection matrix TZa 
of an ideal normal-metal-superconductor interface and the 
reflection and transmission matrices r, r', t, and t' of the 
non-superconducting region. 



are fully quantum mechanical, hence falling outside of the 
scope of quasiclassics. 

We use the scattering matrix approach to calculate the 
differential conductance of an HS junction, and the (zero- 
bias) supercurrent in an SHS junction. We find that there 
is a remarkable difference between these two observables 
in the single-channel limit. For a single-channel half- 
metal-superconductor junction at zero temperature, the 
linear conductance vanishes at the Fermi level. The con- 
ductance becomes appreciable only if the applied voltage 
is comparable to the superconducting gap A or to the 
Thouless energy of the junction, whichever is smaller. 
The Josephson current, on the other hand, proves to be 
nonzero at zero temperature. The origin of this differ- 
ent behavior is that the Josephson effect contains infor- 
mation about the entire excitation spectrum of an SHS 
junction, whereas the linear conductance is a property 
that requires knowledge of excitations at the Fermi-level 
only. 

The remainder of this article is organized as follows. 
In Sec. ini we outline the key elements of the scattering 
approach and its application to HS junctions with a spin- 
active superconductor interface. In Sees. HII] and |TV] we 
then apply the scattering theory to transport through 
an HS junction and to the Josephson effect in an SHS 
junction, respectively. We conclude in Sec. El 



II. SCATTERING APPROACH 

For a scattering description of the triplet proximity 
effect, we consider half-metal-superconductor (HS) junc- 
tions that consist of a half metal "end" , a spin-active in- 
termediate layer, and a superconductor. The intermedi- 



Figure 2: (Color online) Schematic drawing of an SHS junc- 
tion. In the scattering approach, an SHS junction is seen as 
two opposing (composite) HS junctions, with scattering ma- 
trices TZ' and TZ, respectively. In the calculations of Sec. lIV Al 
scattering phase shifts from the central half-metallic part are 
included into TZ' . 



ate layer may be half-metallic, ferromagnetic, or normal 
metallic. 

The central object in the scattering approach is the 
scattering matrix TZ{e) of the HS junction. It relates the 
amplitudes of excitations at energy e > propagating 
towards the superconductor and excitations propagating 
away from the superconductor at the half-metal end of 
the junction, see Fig. [1] If e is below the superconduct- 
ing gap A, all excitations must be reflected at the in- 
terface with the superconductor. This reflection can be 
of normal type (electron-like excitations are reflected as 
electrons, and hole-like excitations are reflected as holes), 
or of Andreev type (electron-like excitations are reflected 
as holes and vice versa). Both reflection types are con- 
tained in the matrix TZ, which is made explicit by the 
decomposition 



(1) 



where r^c and rhh are matrices that describe normal re- 
flection, whereas roh and The describe Andreev reflection. 
All four matrices have dimension iV, the number of prop- 
agating modes at the Fermi level in H. Note that the 
propagating modes in H are not spin degenerate. Be- 
low, we will use the polarization direction of H as the 
spin quantization axis and refer to the electrons with spin 
parallel to the polarization direction of H as "spin up" . 

Knowledge of the scattering matrix TZ is sufficient to 
calculate the conductance of an HS junction, as well as 
the Josephson current in an SHS junction. The zero 
temperature differential conductance of an HS junction 
reads22^ 



G{eV) = ^Tr 



(2) 



(The factor of 2 accounts for the doubling of the current 
by the conversion of an electron into a hole.) An SHS 
junction can be viewed as two HS junctions opposed to 
each other, see Fig. [51 Denoting the scattering matrix 
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corresponding to the second junction as 72.', the Joseph- 
son current reads^^ 

/ = -^^d^ E - Tl'{lOJn)Tl{lOJn)l (3) 

n=0 

where a;„ = {2n+\)'KkBT are the Matsubara frequencies, 
and 5(j) is the phase difference between the two supercon- 
ductors. 

In principle, the exphcit calculation of TZ requires a 
solution of the Bogoliubov-de Gennes equation for the 
full HS junction. Here, we take a different approach^ 
and express TZ in terms of the scattering matrix S of the 
non-superconducting region - that is, the intermediate 
layer and the half metallic region combined - and the 
reflection matrix TZp^ for Andreev reflection off an ideal 
normal-metal-superconductor interface. Using the same 
block structure as in Eq. ([1]), it reads 

^A = a(e)(^_^^^^_,,^^^ Q j, (4) 

where iVs is the number of propagating spin-degenerate 
orbital modes at the Fermi level at the superconductor 
interface and <T2 is the Pauli matrix acting in spin space, 
(j) is the phase of the superconducting order parameter, 
and 

a{e)^e-''''""'°<^/^\ (5) 
The scattering matrix S has the structure 

where S{e) is the scattering matrix describing the 
scattering of electron-like excitations off the non- 
superconducting region. The scattering matrix S{e) can 
be further divided into transmission and reflection blocks, 

where r describes reflection for electrons coming from H, 
r' describes reflection for electrons coming from the su- 
perconductor interface, and t and t' describe transmission 
from and to H. The matrices r and r' have dimension N 
and 2A^Si respectively. Solving for the total scattering 
matrix TZ in terms of 72.A and 5, one then finds 

Toe = r + a^t' a 2r'* (7 2(1 — a^r' a2r'* (J2)~^t , (8a) 

= ie"^at'a2{l ~ a^r'Var-'CTs)"^* , (8b) 

= -ie'"f'at'*a2{l - aVo-ar'Va)"^ , (8c) 

Thh = r* + a^t'*a2r'a2{l - a'r'Var'fTa)^^*. (8d) 

Here we suppressed the energy arguments; the complex 
conjugate matrices in Eq. ([5]) should be taken at en- 
ergy -e. 



In the scattering matrix approach, a necessary condi- 
tion for the superconducting proximity effect is to have a 
nonvanishing rhc. For an HS junction, having a nonzero 
The is not automatic: In the absence of spin-flip scatter- 
ing in the intermediate layer, an electron coming from H 
is Andreev reflected as a spin-down hole. This cannot re- 
enter the half metallic contact; it is reflected from the half 
metal instead, upon which it is Andreev reflected once 
more to return as a spin-up electron. Andreev reflection 
can occur only if the intermediate layer is spin active, 
that is, its scattering matrix is not diagonal in the spin 
up/down basis of the half- metallic contact. Such anoma- 
lous Andreev reflection, in which a spin-up electron com- 
ing from the half-metallic contact is reflected as a spin-up 
hole, is the key to the triplet proximity effect. Examples 
of spin active layers that make this possible are a ferro- 
magnet with a magnetization direction not coUinear with 
the polarization of the half metal, a normal metal with 
strong spin-orbit coupling, or a half-metallic spacer layer 
with a different polarization direction and thin enough 
that there is nonzero transmission of minority electrons 
through evanescent modes. 

In the next two sections we use the scattering theory 
to calculate the conductance of an HS junction and the 
Josephson current in an SHS junction. 



III. HS JUNCTIONS 

A. General considerations 

The scattering matrix TZ{s) obeys particle-holc sym- 
metry, 

n{e) = Ei72(-e)*Ei, (9) 

where Si is the first Pauli matrix acting in electron- hole 
space. For the special case iV = 1, this symmetry, in com- 
bination with the condition that TZ{s) is unitary, leads to 
the condition that either = or Teh = at the Fermi 
level s — 0. As we show in the Appendix, generically 
one has reh(O) = 0, although the possibility roc(O) = 
does occur for certain special choices of the spacer layer. 
The case A'' = 1 is relevant for the case that the con- 
tact to the half metal has only one propagating mode at 
the Fermi level or, alternatively, for the case that there 
is perfect translation symmetry in the transverse direc- 
tion so that different orbital modes do not mix. To the 
best of our knowledge, the observation that Andreev re- 
flection at the Fermi level is absent for single-mode HS 
junctions has not been made before. It presents a quali- 
tative difference compared to FS junctions in which both 
spin directions can propagate. 

In the general theory of Sec. |TT] the spin quantization 
axis is taken to be the polarization direction of the half 
metal. Fixing the spin polarization axis still allows for 
rotations around that axis. For the scattering matrices 
appearing in the theory, such a rotation is represented by 
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the transformation 

1^ e''^'"^/2 J \ e-*''"'='/2 j ' U^J; 

where 5 is the scattering matrix of the non- 
superconducting region, see Eq. the bfock structure 
is that of Eq. ([7]) , and ip is the (azimuthal) angle of the 
rotation. Substituting this transformation into the ex- 
pression ([8]) for TZ, one concludes that such a rotation 
has the same effect on 7?. as a change of the supercon- 
ducting order parameter (/) as 

4>~^(p + ij. (11) 

A consequence of this observation is that, if the inter- 
mediate layer is ferromagnetic or half metallic with a 
polarization along the unit vector 

m — (sin 6' cos ■(/;, sin 6* sin?/), cos 6')"'", (12) 

which makes an angle 9 with the polarization direction of 
the half- metallic contact, 7?. is a function of the difference 
(j) — Tp only. (Here, and in what follows, the polarization 
of the half metal is taken to be along the z axis.) This 
observation, which will be important in our discussion 
of the Josephson effect in SHS junctions below, was first 
made by Braude and Nazarov, using the quasiclassical 
approach.— Here, it appears as a natural consequence of 
the transformation rules of the scattering matrix under 
rotations. 



B. Ballistic HS junction with ferromagnetic spacer 

As a first and simplest application of the theory, 
we consider an HS junction for which the intermediate 
layer is a ferromagnet. The ferromagnet's magnetization 
points along the unit vector given in Eq. (|12p . We take 
the interfaces on both sides of the ferromagnetic spacer 
layer F to be ideal and assume that the electron motion in 



F is ballistic. In that case, different orbital modes decou- 
ple, and one can use an effective single-mode description 
for each orbital mode /i separately. We also assume that 
the thickness of F is short in comparison to the supercon- 
ducting coherence length ^5 — hvp/^ {vp is the Fermi 
velocity), so that the energy-dependence of the scatter- 
ing matrix S can be neglected, and we assume that the 
magnetic flux through F is small in comparison to the 
flux quantum, so that the orbital motion is time-reversal 
symmetric. 

For this system, the calculation of S requires the com- 
position of the 4x4 scattering matrix of the ballistic 
ferromagnetic spacer layer, 

Sf^{^Ij^^, U = e»(''+P™-)/2, (13) 

and the 3x3 scattering matrix Sh of the ideal interface 
between the half-metallic contact and the ferromagnetic 
spacer layer, 



Sh=\1 Q . (14) 



In the above expressions, cr is the vector of Pauli matrices 
(acting in spin space), p — v-^ — vi is the difference of 
the phase shifts of majority and minority electrons in F 
upon propagation through the spacer layer, and 77 = + 
v^. In Eq. ([M]) . /3 is the phase shift spin-down electrons 
experience upon reflection from the half metallic contact. 
The three phases p, 77, and [3 depend on the orbital mode 
p. We have suppressed the mode dependence here, but 
will restore it in the final expression, Eq. p6|) below. The 
block structure of Sp is as in Eq. ([7]). The same is true 
for SH^ where the lower right 2x2 submatrix corresponds 
to the lower right block in Eq. ([7]). 

Combining Eqs. lfT3|) and to calculate 5, and then 
using Eq. ^ to find TZ, we obtain 



7^(£) 



1 + sin^ p siv? ( 



-i/3. 



cos p + i sin p cos 6')^ 



-2i(e/A)e*('^-'^)sin6lsinp e 
I 



2f(e/A)e*(*-''') sinffsinp' 
''•^ (cos p — i sin p cos 0)^ 



(15) 



Substituting Eq. ([S]) for a and summing over all orbital 
modes /i, we conclude that the differential conductance 
of a short ballistic HFS junction is 

G(^) - E¥ (16) 

sir? 6 sir? 

^ A2 ( 1 - sin^ 9 SIT? ) 2 + 4e2 sin^ sin^ ' 



where the summation is over the orbital modes in the HS 
junction. 

This simple result illustrates the two main properties 
of the triplet proximity effect in HS junctions: First, An- 
dreev reflection is possible as soon as there is a spacer 
layer that breaks spin-rotation symmetry around the 
half-metal's polarization direction, provided the elec- 
tron's spin precesses by an angle different from or tt. [In 
Eq. this translates to the requirement that sin 0^0 
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Figure 3: The subgap differential conductance G versus the 
applied voltage V for a ballistic single mode HS quantum 
point contact. The small grey rectangle in the contact repre- 
sents a region with a different magnetization than in the half- 
metallic part. Physically such a region can be present due 
to a misaligned magnetization at the half-metal surface.«i4 In 
our calculations this corresponds to the ferromagnetic spacer 
layer. The curves correspond to different values of the phase 
angles in the ferromagnetic spacer, 9 = 0.8 and p = 0.9 
(dashed curve), 9 = 1.4 and p — 1.2, (dash-dotted curve), 
and 9 = 1.56 and p — 1.53 (dotted curve). 



and sin ^ 0.] And, second, in the absence of or- 
bital mode mixing, G = at the Fermi level, except for 
very special choices of the thickness (proportional to p^) 
and magnetization direction of the spacer layer. In the 
present case, these special choices are angles 9 and 
for which sin^ 9 ~ sin^ = 1. In that case, one finds 
G = {2e'^/h)M, where M is the number of modes with 
sin^ = 1- 

Unlike the quasiclassical approach, the scattering ap- 
proach can also deal with systems in which the num- 
ber of orbital modes is small. The simplest way to il- 
lustrate this is to consider the contribution of one or- 
bital mode; in this case, the result in Eq. and the 
corresponding term in Eq. (jl6p describe a single mode 
ballistic ferromagnetic quantum point contact between 
the half-metal and the superconductor. In Fig. [3] we 
show the differential conductance of such an HS quantum 
point contact for a few representative values of the fer- 
romagnet parameters p and 9. Both aforementioned fea- 
tures are clearly seen: the conductance decreases as sin^ 9 



and sin^ p decreases, and it vanishes at the Fermi energy. 



C. Ballistic HS junction with half-metallic spacer 

If the spacer layer between the half-metallic reservoir 
and the superconductor is not a ferromagnet, but a half 
metal, transmission through the minority channel is via 
evanescent modes, not propagating waves. The scatter- 
ing matrix of the spacer layer, which was given by Eq. 
(fT5)l for the case of a ferromagnetic spacer, now reads 



S 



H' 



(17) 



where 



H' 


















The (mode-dependent) phase shift V]^ and transmission 
coefficient r for minority electrons are functions of the 
wavefunction decay rate q and effective mass of the 
evanescent minority electron wavefunctions, the veloc- 
ity V of the majority electrons, and the thickness d of 
the half-metallic spacer layer, li qd ^ 1, the minority 
electron phase shift i^i becomes independent of the layer 
thickness d, 



ihq/mi 



ihq/mi 



(19) 



whereas the transmission coefficient r oc e"'^''^ and = 
m^vd/h, where is the effective mass of majority elec- 
trons. [The phase shift j3 is the reflection phase for mi- 
nority electron reflection off a half-infinite half metal, see 
Eq. m]) above.] 

With the definitions p — — and rj = + u^, v^e 
then find that the conductance of an HS junction with a 
half-metallic spacer is 



e'^I^^T^ [sinp^ + (1 



(20) 



where we abbreviated 



Bq [sinp^ cose* + (1 - T^)^/^ sin ?7^J^ -f [cos + (1 - r^)^/^cos?7^J^, 
Si = 2^{\^co^^9)(\-T^,)^2{\~T^,fl'^co^{^^,~py){\^CQ^9), 
B2 = 1 + (1 - T^)^/^ cos(r/p - Pp)(l + cos 61) + (1 - r^) cos 61, 

and restored the summation over the orbital modes p. For close to unity, this expression simplifies to the Andreev 
conductance for an HS junction with a ferromagnetic spacer, Eq. (|16p above. For small energies one may neglect the 
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terms proportional to and in the denominator, and we find that G oc e^r. Since the transmission coefficients 
are exponentially small if g^d 3> 1, the conductance is dominated by the transverse mode fx with the lowest q^. 

Similar to the case of ideal transmission, there is a special set of parameters at which the conductance becomes 
large, independent of transmission. This occurs when the coefficient _Bo = in Eq. (|20p . so that the denominator in 
that equation vanishes at e = 0. The condition Bq ^ translates to 

cosp^ = -(1 - T^)^/^ cos 77^, 
sinp^ cos6' = -(1 - T^)^/^sin77^. (21) 

Solutions of Eq. (|2ip satisfy the relation sin^ sin^ d = t^, which generalizes the condition for resonance found for a 
ferromagnetic spacer layer (corresponding to = 1). Since v^^ — (77^ — p^)/2 is a material property if q^d ^ 1, see 
Eq. (fT9|) above, and 77^ are not independent in that limit. For a specific half metallic material and in the limiting 
case q^d ^ 1, the relevant solution of Eq. ([2T|) then becomes {p^ + 77^)/2 = = 7r/2 mod tt and 9 ^ n. Since 
is a function of the thickness d of the spacer layer, not a material property, this condition can always be satisfied for 
special values of d. If a mode satisfies the conditions ([2T|) . its contribution to the conductance is 



2e2 4A2t2 
^ V (22) 

h 4A2T2+£2|[i_cos6i + r^(l + cos6i)]2-4r2}' ^ ' 

At zero energy, one finds perfect Andreev reflection irrespective of t^j. As before, the contribution of a single orbital 
mode in Eqs. (j20p . (|22p describes the differential conductance of a single mode quantum point contact with a misaligned 
half metallic surface layer at the constriction, the analogue of the setup sketched in Fig. [31 



D. Chaotic HS junction 

As the next application of the scattering method, we consider a "chaotic HS junction" , which consists of a half- 
metallic contact, a chaotic quantum dot, a ferromagnetic contact, and a superconductor, all connected in series, see 
Fig. m To illustrate the strengths of the scattering matrix approach, we focus on a situation that is intractable with 
quasiclassical methods: we restrict our discussion to the case that both contacts have one orbital mode only. 

For definiteness, we take the quantum dot to be half metallic, with the same polarization direction as the half- 
metallic contact. Two alternative scenarios, a normal-metal quantum dot and a ferromagnetic quantum dot with a 
magnetization direction parallel to that of the half- metallic contact, will be addressed at the end of this section. In 
all cases we assume that the typical electron path length before exiting from the dot through one of the contacts is 
short compared to the superconducting coherence length. 

The calculation proceeds similar to that of the ballistic junction shown above. For the chaotic HS junction with a 
half-metallic quantum dot, we replace the scattering matrix Sh of Eq. (|14p by 

f-e'^VT^ e^fa+ C)/2^ \ 
Sh = e*(x+«)/2^ e^«VT^ , (23) 
V e'V 

where < t < 1 is the transmission coefficient of the quantum dot and f and x are scattering phases for reflection off 
the quantum dot. As before, /3 is the phase shift minority electrons acquire when they are reflected off the half metal. 
The special case r = 1 simplifies to the ballistic HS point contact we considered previously. For general r, one finds 

^, , 2e2 lee^AV^sinVsin^^' 

= - w) ' ^''^ 

where 

D = I A2(1 - t)1/2 [4 (jQg g gjj^(^ _ ^) sm(2p) + cos(/3 ^ ^) {l - 2 cos(26l) sin^ p + 3 cos(2p))] 

+ 2(r - 2) (2e2 + A^ sin^sin^ 9 - A^) }' + We^ (A^ - e^) t\ (25) 

I 



For the special point sin^ = sin^ p = 1 one has G'(O) = 2e^//i, independent of t and the scattering phases 
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Figure 4: (Color online) Chaotic HS junction, consisting of a 
half- metallic contact (left), a half-metallic quantum dot (cen- 
ter), a ferromagnetic spacer layer, and a superconducting con- 
tact (right). 



(3, X, and ^. The origin of this remarkable result is that, 
for sin^ 6 = sin^ p = 1, the ferromagnet-superconductor 
interface not only provides perfect Andreev reflection be- 
tween spin-up electrons and spin-up holes, but, moreover, 
after two subsequent Andreev reflections the net phase 
shift is — = 1 at the Fermi energy. Hence, combining 



the interface reflection matrix (jlSp with the scattering 
matrix of the quantum dot, the conductance at vanish- 
ing voltage is found to be 



2e2 



G(0)=rV[l-(l~r)]^ 



1, 



(26) 



independent of the dot's transmission coefficient t. 
This is to be contrasted to the corresponding formula 
formula— 



(27) 



for the linear response conductance (per spin) of a su- 
perconductor in contact with a single mode quantum dot 
through a normal-metal contact (and without magnetic 
a field). The difference arises from the fact that, in the 
latter case, the phase shift upon two Andreev reflections 
is = — 1 at the Fermi energy. 

For a chaotic quantum dot, the transmission coefficient 
T and the scattering phases < ^, x < 47r are random 
quantities with the statistical distributioni^ 



1 



327r2 



-1/2 



(28) 



As is standard in the statistical approach to quantum 
transport, the statistical ensemble is obtained by means 
of small variations of the dot's shape or of the Fermi 
energy. In an experiment, both types of variations can be 
achieved by changing the voltage of nearby metal gates. 
With the help of the distribution ([^5]) we can calculate 
the average Andreev conductance (G) for an ensemble 
of quantum dots. The phase shift p, the angle 9 and 
the reflection phase /3 are not averaged over, since they 
are properties of the ferromagnetic contact and the half 
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Figure 5: The ensemble averaged subgap differential conduc- 
tance (G) versus the applied voltage V for different phase 
angles 9 and p describing the ferromagnetic contact. The val- 
ues of 6 and p are 0.8 and 0.9 (dashed curve), 1.4 and 1.2, 
(dash-dotted curve), and 1.56 and 1.53 (dotted curve), respec- 
tively. The solid curve shows the special case corresponding 
to 6* = p = 7r/2- 



metal interface, not of the chaotic quantum dot. The 
average can be performed analytically in the special case 
sin^ p — sin^ & — ^, for which we find 



with 



where 



2e2 

^'^^"Xl^ 4e(£2-A2)i/2 



S = ^ zbx^artanh— , 



S 



2£(e±(e2_ A2)i/2) 



(29) 



(30) 



(31) 



At the Fermi level (G) = 2e^/h, in agreement with the 
discussion following Eq. (I25|) . For general values of 9 and 
p no closed-form expression for (G) could be obtained. 
The result of a numerical evaluation of the ensemble av- 
erage (G(e)) is shown in Fig. [5] for a few representative 
values of and p. For generic 9 and p, the ensemble- 
averaged conductance vanishes at the Fermi level e — 0. 
The quadratic dependence G ~ for e ^ changes to 
a linear increase for relatively small voltages. The con- 
ductance reaches a maximum at a voltage eV below the 
superconducting gap A for 9, p sufficiently away from 
sin 9 — 0, sin p = 0. With 9, p approaching sin 9 — 0, 
sin p = 0, the position of the maximum moves towards 
the superconducting gap A. For 9, p close to sin9 — 0, 
sin p ~ 0, the conductance is an increasing function of 
the voltage in the full subgap regime. 

For the case of a normal- metal quantum dot, the con- 
ductance G(e) is given by Eq. ([M)) . but with the replace- 
ment (3 — ^ ^ arctan[T sin(/9 + x)/((2 — ''') cos(/3 + x) + 
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2(1 — tY^"^)]- For the case of a ferromagnetic quantum 
dot with a magnetization direction along that of the po- 
larization of the half-metallic contact, G{e) is given by 
Eq. (|24p . but with /? — ^ replaced by a phase shift /3' that 
is statistically independent of x, ^, and r. In both cases, 
the qualitative dependence of (G) on the parameters e, 
6, and p is the same as in the case of a half-metallic 
quantum dot discussed above. 

IV. SHS JUNCTIONS 

We now contrast the transport current through an HS 
junction to the supercurrent through an SHS junction. 
As in the previous section, we consider the effect of a 
thin ferromagnetic layer between each superconductor 
and the adjacent half metal. (We do not consider the 
case of a thin half- metallic spacer layer in this section.) 
While, at zero temperature, the zero-bias conductance of 
a single single-channel HS junction vanishes (except at 
special choices of the parameters) , the zero temperature 
Josephson current / is not zero. The reason is that, in 
contrast to the linear response conductance G, I is not 
a Fermi level property. Instead, it is determined by the 
full excitation spectrum. 

In order to apply the theory of the previous sections, 
we consider the SHS junction as two opposing HS junc- 
tions, see Fig. [2] We refer to the opposing HS junction 
as S'H. Both junctions have intermediate ferromagnetic 
layers, which are denoted by F and F'. The two ferro- 
magnets can have different magnetizations, parameter- 
ized by polar angles 9, tp and 0' , ip' , respectively. The 
superconductors S and S' are assumed to have equal su- 
perconducting gaps A, but the phases (j), 0' of the order 
parameters can differ. 

Before turning to applications of our scattering theory, 
it is worthwhile to summarize some general considera- 
tions. Because of the transformation property pT|) . the 
Josephson current / can depend on the superconducting 
phases cf) and (p' and on the azimuthal angles V' and ip' 
through the single combination 

^ = - 0' - (^ - V') (32) 

only. This observation was made previously in the con- 
text of the quasiclassical approach. ^Siiiii^ 

Under the operation of time reversal, the phases of the 
superconductors and the (position dependent) magneti- 
zation direction m transform a,s cj) ^ — (/),m^ — m. The 
supercurrent of the time reversed system is the opposite 
of the original, that is, 

/(0-0',m) = -/(<^'-0,-m). (33) 

The supercurrent is invariant under a position indepen- 
dent rotation of the magnetization. This, together with 
Eq. (I32D results in /(0) = -I{-4>). 

For phase angles not close to the special point sin^ 9 = 
sin^ 6*' = sin^ p = sin^ p' = 1, the Andreev reflection 



probability at the SH interfaces is significantly smaller 
than unity [see Eq. pS)) above]. As a consequence, the 
(/)-dependence of the supercurrent is nearly sinusoidal in 
this case. The detailed calculations of the next section 
show, however, that close to the special values of the 
phase angles the ^-dependence becomes non-sinusoidal. 

As an illustration of our scattering theory, we now con- 
sider the ballistic and chaotic junctions addressed in the 
previous section. Our work on the Josephson effect in 
ballistic junctions complements that of Galaktionov et 
al., who used a Green function approachji^ 

A. Ballistic SHS junction 

For the ballistic SHS junction different orbital modes 
are not mixed, so that the scattering problem is effec- 
tively one-dimensional. As before, we denote the differ- 
ence of the (mode-dependent) phase shifts of majority 
and minority electrons transmitted through F by p, see 
Eq. (|13p: The corresponding quantity for F' is denoted 
by p' . We suppress the mode index p, except in the fi- 
nal expressions. For the calculation of the supercurrent, 
it is necessary that phase shifts accumulated inside the 
half metal are included into the determinant in Eq. ([3]). 
For an orbital mode p these phase shifts depend on the 
length L of the half-metallic segment and on the longi- 
tudinal component kf^(s) = fc^(O) -t- e/{hv^j) of the wave 
vector for that mode, where Vn is the group velocity of 
the mode at /c^(0). In order to include this into Eq. ([3]) 
we take the scattering matrix TV to include the scattering 
phase shifts accumulated inside the half metal, 

\ f^~ik^{~e)L j '"^ y g-ik^(~<^)L j ' 

(34) 

where TV is the reflection matrix for the S'H junction 
without the scattering phases from the half metal. This 
matrix is given in Eq. (|15p of the previous section, but 
with 9, tp, (p, and p replaced by 9' , tp' , cp' , and p', respec- 
tively. 

Since there is a probability of normal reflection at each 
end of the SHS junction, for a given orbital mode, the 
contribution to the supercurrent contains terms that os- 
cillate with the length L of the junction. For the total su- 
percurrent, obtained by summing the contributions from 
different orbital modes, however, this results only in a 
small correction, provided that /c^(0)-L 3> 1, since in this 
case, the sum of the oscillating contributions averages 
out. Below, we calculate the non-oscillating contribution 
to the Josephson current for a given orbital mode, and 
restrict our discussion to the limiting cases of a "short 
junction" [L <C Cs) and a "long junction" (L » ^5). 
(In both cases, we assume that the ferromagnetic spacer 
layers are thin in comparison to the superconducting co- 
herence length ^5. The same assumption was made in 
the previous section.) 



-0.6 



-1 



Figure 6: (Color online) The contribution of a single trans- 
verse mode to the zero temperature supercurrent / of a short 
SHS junction, as a function of (j), for ferromagnetic phase an- 
gles e = e' = p = p' = Tx/2 (solid), e = e' ^ p = p' = Tx/A 

(dot-dash), and 6 = 9'= 7r/2, p = p' = tv/A (dashed). The 
supercurrent is shown in units of /short = eA//i. 



For a short junction, one may neglect the energy- 
dependence of the wavenumbcr fc^(e) in the half metal. 
A closed-form expression valid for arbitrary temperatures 
could be obtained for the special case sin^ 9 = sin^ 6' = 1 
for a mode fx with sin^ = sin^ pji = 1 only. The con- 
tribution to the supercurrent of such a mode is 



eA ip + s^-K , / A . ip + s^Ti 

In = cos ^— tanh sm ^— 

2h 2 \ 2kBT 2 



where is defined through the relation 
(-1)"^ = sinp^^sinp^. 



(35) 



(36) 



The TT shift in the current-phase relationship associated 
with originates in the properties of the interface re- 
flection matrix (jlSp : for this matrix, the transformation 
p — > p -|- TT is equivalent to (p ^ (j) + n. 

In the limit of high temperatures fc^T 3> A, one can 
find a closed- form expression for arbitrary values of 0, 9' , 
p, and p' . Upon summation over all orbital modes, one 
has 



i-i- 



T 



sm cp sm sm p^ sm 9 sm ( 



(37) 



Note that although the angles p^, p^ are mode depen- 
dent, for sufhciently thin spacer layers the mode depen- 
dence is weak enough that all modes contribute to the 
total Josephson current with the same sign. The super- 
current is reduced once the thickness of the spacer layers 
is large enough that p^, p^ ^ 1. 

A numerical evaluation of the contributions to the zero- 
temperature supercurrent is shown in Fig. [S] for a few 
choices of the_angles 9, 9', p, and p'. Although the dis- 
continuity at (j) = sir is smeared for generic values of the 
phase angles, the order of magnitude of the supercur- 
rent is the same as at the special point sin^ 9 = sin^ 9' = 



Figure 7; (Color online) The contribution of a single trans- 
verse mode to the non-oscillating component of the zero tem- 
perature supercurrent / of a long SHS junction, as a function 
of (j), for ferromagnetic phase angles 6 = 6' = p = p' = n /3 
(solid) and 6 = 9' = p = p' = Tr/4: (dashed). The current 
is shown in units of /long = eft^n^/vrL'^A^, 
mode-dependent longitudinal velocity. 



sin^ p 



sin^ p' 



1. This is in contrast to the Fermi- 
level conductance of an HS junction, which was zero for 
generic phase angles and finite at the special point. As 
discussed above, the reason why the supercurrent has a 
different behavior is that it is not a Fermi-level property 
but, instead, depends on the entire excitation spectrum. 
For energies far away from the Fermi level, the Andreev 
conductance is not qualitatively different at the special 
point and elsewhere, see Fig. [3l 

For a long SHS junction (but still with ferromagnetic 
spacer layers that are much thinner than ^5), again a 
compact expression at arbitrary temperatures could be 
obtained for the special case sin^ 9 = sin^ 9' = 1, for the 
contribution of a mode p with sin^ pfj, = sin^ p'^ = 1 
only. In this case one finds 



sin(^ + s^tt) 



cosh{2ujnL/hv^) - cos(0 -I- s^tt) 



where was defined in Eq. (j36p . At zero temperature 
the summation can be replaced by an integration and one 
has 

ZttL 

In the limit of high temperatures, T 3> hVf^/L, only the 
term with n = contributes, so that 

= -liksTe-^^''^^''^^^- sin(0 + s^t:). (40) 

The special point sin^ 9 — sin^ 9' = sin^ p — sin^ p' — ^ 
is singular, however, and the supercurrent contributions 
have a qualitatively different dependence on temperature 
for generic 9, 9', p, and p'. In the high-temperature 
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Figure 8: (Color online) Superconductor-half-metal- 
quantum-dot -superconductor junction. In the calculation, 
scattering from the half-metal quantum dot is included in 
the scattering matrix TZ. 



regime hv^/ L ^ fc^T ^ A, one finds 



167r2/c3r3 



E 



sin sin p'^ sin 9 sin 0' sin ( 







o 



-0.25 











TT 



Figure 9: The ensemble averaged Josephson current (/) as the 
function of (j> for different phase angles 9, p, 6' , p' describing 
the ferromagnetic spacer layers. The values of 6 and p are 
e ^ p = Tr/2 (solid curve), 0.9 n/2 and 0.99 7r/2 (dashed 
curve), and 0.5 7r/2 and 0.6 tt/2 (dotted curve), respectively. 
The values for the second contact are 0' — 1.05 n/2 and p' — 
0.95 7r/2. The supercurrent is shown in units of /short ~ eA/h. 



'2TTLkBT/hv^ 



(l — sin^ p^ sin^ I 



1 



sin^ p;, sin^ 0') 



(41) the quantum dot. Since (r^/^) — 1/2, at zero tempera- 
ture, the ensemble-averaged supercurrent is 



This result is a factor ~ (fc^T/A)^ <C 1 smaller (per or- 
bital mode) than the contribution for the special choice 
of the angles 9, 9', p, and p' in Eq. (j40|) . Whereas the 
supercurrent of a short Josephson junction depends on 
the full subgap excitation spectrum of the junction,^- 
the supercurrent in the long junction limit is determined 
by the junction's excitation spectrum up to the Thouless 
energy hvp/L only.^^ In this range of the spectrum, the 
absence of Andreev reflection at the Fermi energy still 
strongly affects the magnitude of the supercurrent. For 
temperatures below the Thouless energy hvp / L the sup- 
pression factor with which is reduced in comparison to 
the special case of Eq. ([55)) saturates around (Hvp/LA)'^. 
No closed-form expressions for at arbitrary tempera- 
tures could be obtained. Figure [7l shows versus at 
zero temperature, for two choices of the parameters 9, 9' , 
p, and p' . 



(I) 



eA 

I 

4n 



< (p + STT < 2tt, 



(42) 



where s was defined below Eq. ([35]) . No closed- form ex- 
pressions could be obtained for generic values of 9, 9', p, 
and p' . A numerical evaluation of the ensemble-averaged 
supercurrent is shown in Fig. [9l For phase angles close 
to the special case discussed above, the supercurrent in 
a short chaotic SHS junction follows Eq. except 
near the discontinuity (f> = stt, which is smoothed out 
away from the special point sin^ 9 = sin^ 9' — sin^ p — 
sin^ p' = 1. For phase angles not close to the special 
point, the ^-dependence of the ensemble- averaged super- 
current is nearly sinusoidal, confirming the general ob- 
servations made in the beginning of this section. 



CONCLUSION 



B. Chaotic SHS junction 

For a chaotic SHS junction, we include the quantum 
dot into TZ, and take TZ' — IZ' to be the scattering matrix 
of a junction without quantum dot, see Fig. [S] For the 
chaotic SHS junction we only consider the limit that the 
superconducting coherence length is much longer that 
the typical electron path length in the dot before exiting 
through one of the contacts. 

In the special case sin^ 9 — sir? 9' = sin^ p — sin^ p' — 
1 the expression for the supercurrent is the same as 
for a ballistic SHS junction, but with the replacement 
A At^/^, where r is the transmission coefficient of 



For the conventional proximity effect, the possibility 
of Andreev reflection of electrons at the Fermi level gives 
a nonzero linear conductance through a normal-metal- 
superconductor interface. In this article, we found that 
the situation is more delicate for the triplet proximity 
effect in half-metal-superconductor (HS) junctions. In 
the case that there is only one conducting channel at 
the HS interface, or that different orbital channels at the 
HS interface decouple, we found that Andreev reflection 
processes can be present only away from the Fermi level 
(except for special choices of the interface parameters). 
While this result, which is independent of the nature of 
the spin active spacer layer in the HS junction, leads to 
a vanishing linear conductance, it allows for a nonzero 
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Josephson current through an effectively single-channel 
SHS junction. We have illustrated this statement on sys- 
tems both in the quasiclassical and in the fully quantum 
mechanical regimes. In our calculations we have mainly 
concentrated on the case of ferromagnetic spin active in- 
termediate layers. 

First, we have calculated the zero temperature differ- 
ential Andreev conductance at finite bias for short HS 
junctions. This is the observable in which the present 
Andreev reflection processes manifest themselves in the 
most direct way. Using the scattering matrix approach, 
we calculated the dependence of the Andreev conduc- 
tance on the phase angles of the spacer for all subgap 
voltages. Our result in Eq. ([M]) can be used to describe 
the conductance of a system with an arbitrary single 
channel structure in the half metal, provided that its 
normal state scattering matrix is known. As an appli- 
cation, we considered the case that the structure is a 
chaotic quantum dot and we calculated the ensemble av- 
eraged conductance from the known distribution of the 
dot scattering matrices. In addition to the calculation of 
the differential conductance for systems with ferromag- 
netic spacer layer, we also studied ballistic systems where 
the spacer is a thin half metallic layer. 

Second, we calculated the zero bias Josephson current 
through SHS junctions. We have confirmed the observa- 
tion, reported in earlier works jiSiiiiii^ that the Josephson 
current depends on the superconducting phase through 
the single variable = (f> — (/)' — {iJj — ip') only, which 
is the difference of the superconductor phase difference 
of two superconducting reservoirs and the azimuthal an- 
gle differences of the magnetization direction of the two 
ferromagnetic spacer layers in the SHS junction. In the 
framework of the scattering matrix approach, this obser- 
vation follows directly from the fact that the phase of the 
superconductor and the azimuthal angle of the ferromag- 
netic spacer at an HS interface enter in identical ways 
in the calculation of the Andreev reflection amplitude. 
Further symmetry considerations showed that the super- 
current is an odd function of the variable 4>. Similarly 
to earlier worksiSiiiiiiiii^, we also find that for symmetric 
ferromagnetic spacers, ip = ip' , p ^ p' , 9 ~ 9', the sign 
of the current is the opposite to the case of conventional 
SNS junctions [see Figs. I6I7I5] . Consequently, the equi- 
librium phase difference corresponds to (p — cj)' = tt , i.e., 
a TT-junction behavior is realized. For independent con- 
figurations in F and F', the equilibrium phase difference 
varies continuously as the function of interface parame- 
ters. 

It is worthwhile to compare our results for the Joseph- 
son current in single channel SHS systems to the result 
for a single channel SNS systems. In the latter case, at 
zero temperature and in the absence of magnetic field, 
for a perfectly transparent normal region, the (per spin) 
Josephson current is given hy^ I = {eA/2h) sin((^/2), for 
short junctions and^^ / = ewi?'/'/27rL for long junctions, 
where < tt. We found [see Eqs.dSlD and ^] that 
in the case of single channel SHS systems, in the special 




Figure 10: Sketch of a possible experimental setup for testing 
the vanishing Andreev reflection at the Fermi level: a single 
channel quantum point contact to an FS junction. The arrow 
in the quantum point contact indicates that the point contact 
transmits only one spin direction. 



point sin^ 9 — sin^ 9' — sin^ p ~ sin^ p' = 1, the current- 
phase relation is identical, apart from the phase shifts 
due to the azimuthal angles and s. Away from the spe- 
cial point, the current-phase relation becomes sinusoidal, 
similar (apart from the phase shifts) to the case of a 
normal region with low transparency. By adjusting the 
interface parameters, the single mode triplet Josephson 
current interpolates between the result for the conven- 
tional Josephson current through an ideal single mode 
channel and through a tunnel barrier. The key prop- 
erty that distinguishes the current phase relation in the 
triplet Josephson effect through single mode structures 
from the conventional Josephson effect is the magneti- 
zation dependent phase shift. This is a feature that is 
common between the fully quantum mechanical single 
channel limit and the multi mode case corresponding to 
the quasi-classics. 

We end by relating our results about HS junctions to 
a possible experiment. One experimental setup could be 
the HS quantum point contact sketched in Fig.[3l Such a 
setup is somewhat subtle, as it relies on the presence of a 
surface magnetization in the point contact. The general- 
ity of our proof in the Appendix suggests, however, that 
the main features of the single channel HS conductance, 
i.e., G = at Fermi level and G 7^ for < el/ < A could 
be tested in an experimentally more robust arrangement. 
Such a setup could be a single channel point contact to 
an FS junction, as sketched in Fig.[TOl It is not necessary 
to have the system in the short junction limit, and there 
can be arbitrary number of modes at the ferromagnet- 
superconducting interface. The only important detail is 
that the junction ends in a single mode point contact 
through which only one spin direction can be transmit- 
ted. This can be achieved using a half metallic electrode 
or with a spin filtering quantum point contact^. 
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APPENDIX: ABSENCE OF ANDREEV 
REFLECTION FOR SINGLE-MODE HS 
JUNCTIONS 

In this appendix we prove that, generically, the An- 
dreev reflection amplitude rhe(O) = for a junction with 
= 1 orbital modes in the half metallic side. The 
number of modes on the superconducting side can be 
arbitrary. The starting point of the proof is the sin- 
gular value decomposition of the scattering matrix S of 
the non-superconducting region between the half-metallic 
and superconducting reservoirs, 

-(o:)(fS)(o ;.)■ <«' 

Here, V and V' are N x N unitary matrices, W and 
W are unitary matrices of dimension 2A^5, Ag being 
the number of orbital channels at the normal-metal- 



superconductor interface, T is an 2A^s ^ N matrix with 
fki=SkiVn, k = l,...,2Ns, l^l,...,N, (44) 
with Ti the ^th transmission eigenvalue, ^ = 1, . . . , A, and 

R=\/1n- fTf, R'^^ t2Ns - TfT. (45) 

Substituting the decomposition in Eq. ([5c)) . and as- 
suming det(l27Vs + r' ij2r'* 1J2) 7^ 0, one finds 

rhe(O) = -^e-''^V*f^ {z^ - R'Z*R'^ ^ TV' (46) 

with Z = W'*a2W. If A^ = 1, the amphtude rho(O) 
is proportional to the 11 element of the inverse in (|46p . 
Using the general result A^^ = (dot A)^^adj(^) for the 
matrix inverse, we find that this element is proportional 
to the determinant of an antisymmetric matrix of dimen- 
sion 2A^s — 1, and is therefore zero. The case rhc(O) 7^ is 
possible if det(l2Ars -t- r'a2r'*(T2) = 0, that is, if the sys- 
tem has an Andreev bound state at e = that is not 
coupled to the mode in the half metal. For the ballistic 
HS system in Sec. IhTbI 

det(l2JVs + r'(j2r'*(J2) = 1 - sin^ psin^ 6*, (47) 

resulting in sin^ p = src? = 1 to be the only points 
where rhc(O) can be nonzero. 
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